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摘 要 :直升机 航 炮 在 发 射 时 形成 的 强烈 冲击 波 不 仅 对 周边 结构 构成 损伤 ,产生 的 炮 振 环境 还 会 导 
致 附近 设备 的 失效 和 损坏 , 航 炮 的 机 弹 相 容 设 计 一 直 是 武装 直升机 研制 过 程 中 的 难点 问题 。 本 研 
究 基 于 经 典 的 爆炸 力学 理论 和 某 航 炮 地 面 靶 试 实测 结果 ,对 炮 口 冲击 波形 成 的 压力 场 进 行 了 数值 
拟 合 , 同 时 考虑 冲击 波 正 压 时 间 的 差异 性 ,分 别 建立 了 超 压 分 布 的 空间 场 与 时 间 脉 冲 场 。 根 据 模拟 
台 架 与 航 炮 的 结构 方案 ,在 有 限 元 软件 中 对 炮 口 冲击 波形 成 的 超 压 扩散 过 程 进 行 了 模拟 ,然后 利用 
非 线性 肯 态 分 析 方 法 ,对 直升机 模拟 台 架 在 超 压 作用 下 的 冲击 过 程 进行 了 仿真 ,分 析 了 人 台 架 蒙 皮 、 
骨架 以 及 连接 结构 的 动态 响应 变化 规律 ,可 以 为 地 面 台 架 的 尺寸 设计 和 结构 优化 提供 依据 。 
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Numerical simulation of helicopter simulation platform 
under the action of aerial gun shock wave 
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;Abstract :The strong shock wave generated by the helicopter gun during launch not only damages the sur- 
rounding structure ,but also causes the failure and damage of nearby equipment ,the compatibility design of 
aerial gun with projectile is always a difficult problem in the development of armed helicopter. Based on 
the classical theory of explosion mechanics and the measured results of an aerial gun ground target, the 
pressure field generated by the muzzle shock wave is simulated numerically ,and the difference of the time 
of the positive pressure of the shock wave is considered ,the space field and time pulse field of overpressure 
distribution are established respectively. According to the structure scheme of the simulation platform and 
the aerial gun ,the overpressure diffusion process of the muzzle shock wave is simulated in the finite ele- 
ment software ,and then the nonlinear transient analysis method is used, the impact process of helicopter 
platform under overpressure is simulated ,and the dynamic response of skin ,frame and connecting structure 
is analyzed , which can provide reference for the design and structure optimization of ground platform. 
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航 炮 作为 航空 器 机 载 的 重要 武器 ,在 近 距 作战 、 
对 地 打击 等 方面 发 挥 着 重要 作用 ,已 经 成 为 现代 武 
装 直升机 的 必 备 装备 之 一 ") 。 航 炮 发 射 形成 的 强烈 
冲击 波 、 温 度 场 和 噪声 波 等 往往 引起 周边 结构 的 站 
劳 破坏 ,其 至 导致 机 载 设备 的 失效 和 损坏 '””。 为 此 ， 
航 炮 的 机 弹 相 容 问 题 已 经 成 为 飞行 器 设计 中 的 重大 
风险 项 之 一 。 关于 航 炮 设计 ,国外 的 研究 相对 较 早 ， 
已 形成 了 完善 的 理论 体系 ,国内 的 相关 研究 起 步 
较 晚 ,早期 主要 采用 数值 计算 和 试验 的 方法 ,对 膛 口 
冲击 波 的 传播 规律 进行 了 研究 。 李 鸿 志 等 采用 高 
精度 Godunov 格式 ,对 膛 口 冲击 波 在 远 场 的 传播 过 
程 进行 了 数值 模拟 ,得 到 了 动 球 心 球形 波 特征 的 冲 
击 波 结果 ,给 出 了 超 压 在 传播 过 程 中 呈 指 数 衰减 的 
规律 。 郭 则 庆 等 “基于 ALE 方程 的 有 限 元 体积 法 ， 
对 瑟 计 算 了 有 无 初始 流 场 的 复杂 形状 弹丸 在 膛 口 的 
流 场 发展 过 程 , 结 果 表 明 初 始 流 场 是 火药 燃气 流 场 
形 雹 冠状 冲击 波 的 必要 条 件 ;随后 ,又 对 高 马赫 数 下 
顶 旨 内 埋 航 炮 流 场 进 行 了 计算 ,指出 高 速 飞行 时 及 
HE 佬 形成 由 火药 燃气 冲击 波 .分离 激 波 、 滑 移 面 等 波 
系 构成 的 流 场 结构 ,分 离 激 波 尺寸 和 冲击 波 超 压 峰 
全 到 与 马赫 数 有 关 !" 。 朱 冠 南 等 ”分 别 采 用 高 、 低 
虑 胎 装 置 , 在 实验 室 进行 了 步枪 膛 口 冲击 波 的 压力 
测 焉 试验 , 测 得 不 同 压力 下 膀 外 多 点 处 的 压力 曲线 ， 
通过 数值 拟 合 , 给 出 了 指数 衰减 的 近似 公式 ,并 对 比 
分 煌 了 环境 压力 对 冲击 波 强度 的 影响 。 
> 泛 随 着 计算 机 和 有 限 元 技术 的 发 展 ,基于 三 维 流 
场 腊 析 的 复杂 工程 结构 冲击 波 结构 响应 分 析 与 优化 
成 海 当 前 的 趋势 ,被 广泛 用 于 制 退 器 ”、 炮 车 、 
装甲 车 、 舰 载 航 炮 "等 相关 结构 "的 设计 。 
但 对 于 航 炮 ,特别 是 直升机 航 炮 设计 ,工程 上 多 依赖 
于 试验 ,缺乏 有 效 的 分 析 手 段 。 

为 了 对 某 航 炮 装机 冲击 波 影响 进行 分 析 , 本 研 
究 首先 基于 文献 中 的 理论 和 经 验 公 式 , 对 实测 炮 口 
冲击 波 的 空间 压力 场 进行 数值 拟 合 , 并 采用 非 线性 
瞬 态 分 析 方 法 ,对 某 直 升 机 模拟 试验 台 架 在 冲击 波 
作用 下 的 结构 响应 进行 数值 仿真 ,基于 分 析 模型 , 开 
展 了 有 关 参 数 的 影响 分 析 , 用 于 台 架 的 尺寸 设计 与 
优化 。 


1 炮 口 冲击 波 分 布 


1.1 超 压 分 布 
国内 学 者 通过 数值 方法 对 膛 口 气 场 的 形成 、 流 
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动 过 程 进行 了 模拟 ,基本 确定 了 炮 口 冲击 波 的 远 场 传 
播 规律 ” ,并 通过 相关 的 枪击 和 炮击 试验 “进行 了 
验证 ,给 出 了 膛 口 冲击 波 超 压 在 不 同 角度 射线 上 呈 震 
函数 衰减 的 规律 ,参考 上 述 文献 , 寡 函 数 形式 如 下 。 
AP = 天 r” (1) 
其 中 :AP 为 超 压 值 (kPa) ;kK 为 起 始 冲击 波 超 压 值 ; 
as 为 超 压 值 衰减 率 ;r 为 测 点 到 炮 口 直线 距离 与 炮 
口 直 径 的 比值 。 
1.2 正 压 时 长 
炮 口 冲击 波 传 播 到 空间 任意 点 的 压力 场 强度 ， 
除了 扩散 角度 的 衰减 外 , 男 一 个 关键 参数 为 正 压 的 


作用 时 间 , 由 爆炸 相似 律 ” 可知, 炸药 爆炸 空中 冲 
击 波 的 正 压 时 长 满足 


t d 
=g| 3 二 2 
志 -a( 浙 | 2) 
式 中 :4, 为 冲击 波 正 夺 作用 时 间 (ms) ;m 为 TNT 装 药 量 
(kg) ;d 为 雹 口 距离 (m) :定义 ;= 为 比例 虐 离 。 
文献 [18] 根 据 实测 的 试验 数据 对 近 地 面 空中 爆 
炸 冲击 波 正 压 时 间 进 行 了 数值 拟 合 , 并 与 其 他 文献 中 
结果 进行 了 对 比 验证 ,给 出 了 更 为 精确 的 关系 式 为 


二 a es 
= -0.74+0.97d-0.04¢4,1<d<10 (3) 


3 
Vm 


2 数值 模型 


2.1 动力 学 模型 


以 某 直升机 模拟 台 架 结构 为 分 析 对 象 , 台 架 由 
骨架 、 蛇 皮 及 设备 模拟 件 组 成 。 其 中 骨架 为 Q235 钢 
制 的 标准 方 管 , 蒙 皮 分 别 采 用 钢板 和 铝板 。 连 接 均 采 
用 焊接 形式 ,设备 模拟 件 为 木 制 材料 ,通过 4 个 螺栓 
安装 在 侣 架 前 部 平台 。 人 台 架 的 结构 模型 如 图 1 所 示 。 


\ 7 
/Nf 

Z 
Ja 
1 一 骨架 ; 2 一 设备 模拟 件 ，3 一 钢 蒙 皮 ; 4 一 铝 蒙 皮 


图 1 直升机 模拟 试验 台 架 结构 方案 


Fig.1 The scheme of test platform for simulating 


helicopter structure 


对 上 述 结构 进行 动力 学 建 模 。 骨 架 采 用 梁 单 
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染 采 用 党 散 。 在 离 炮 口中 心 不 同 射 角 和 距离 处 布置 压力 传 感 


元 , 莹 上 色 和 设备 模拟 件 采用 2D 充 单 元 ,模型 中 单元 
尺寸 为 10 mm。 材料 包括 Q235 钢 和 2A12 铝板 , 考 
虑 结构 大 变形 ,在 模型 中 创建 金属 材料 的 应 力 -应 变 
曲线 ,分 析 时 采用 非 线 性 材料 本 构 。 


2.2 材料 非 线 性 


金属 材料 在 大 变形 过 程 中 可 以 分 为 弹性 段 和 塑 
性 段 , 其 中 , 弹性 段 内 , 应力- 应变 存 在 一 一 对 应 关 
系 ,由 明 克 定律 可 知 


ry NO +0 (4) 
A=vE/(l +v) (1 -2v) (5) 
G=E/2(1 +v) (6) 


其 中 :2 为 拉 梅 系数 ; 6 为 剪 切 模 量 ;E 为 弹性 模 量 ; 
2 为 泊 松 比 。 
69 当 进入 塑性 变形 阶段 ,载荷 印 去 存在 残余 变形 ， 
此 时 ,应 力 - 应 变 之 间 不 存在 唯一 的 对 应 关系 ,曲线 
如 图 2 所 示 。 
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Fig.2 Stress-strain curves of elastoplastic materials 


对 于 不 可 压缩 材料 ,应 变 a 和 真实 应 力 ou 


us =ous(1+e) (7) 
趟 中 款 。 为 工程 应 力 。 

对 于 各 向 同性 的 金属 材料 ,本 研究 中 采用 Von 
Mises 届 服 准则 、 关 联 流动 准则 和 各 向 同性 硬化 定律 
对 材料 本 构 进行 描述 。 

根据 实际 的 试验 边界 ,约束 台 架 后 端面 4 个 接 
头 的 所 有 平 动 自由 度 , 建立 的 动力 学 模型 如 图 3 
所 示 。 


2.3 压力 场 拟 合 


某 航 炮 地 面 靶 试 试验 时 ,采用 地 面 刚性 基础 对 
炮 身 进行 固定 , 炮 口 为 自由 端 ,便于 冲击 波 自由 扩 


器 ,实测 航 炮 发 射 时 的 冲击 波 超 压 峰 值 , 测 点 位 置 如 
图 4 所 示 。 


图 3 直升机 模拟 试验 台 架 动力 学 模型 


Fig.3 Dynamics model of helicopter simulation test platform 


图 4 地 面 丢 试 超 压 测 点 布置 

Fig.4 Location of overpressure sensors in ground test 

对 地 面 靶 试 试验 的 冲击 波 实测 结果 进行 处 理 ， 
得 到 不 同 角度 和 不 同 炮 口 距离 处 测 点 的 超 压 值 。 超 
压 值 与 炮 口 距离 比 和 扩散 角度 的 关系 如 图 5 所 示 。 


0 30 60 90 120 150 
扩散 角度 /9) 


图 5 炮 口 冲击 波 超 压 值 与 距离 比 和 角度 的 对 比 


Fig.3 The overpressure of muzzle shock wave varies 


with the distance ratio and angle 


630 
采用 式 (1) 中 表达 式 , 对 上 图 中 航 炮 发 射 冲击 
波 的 实测 压力 场 进行 最 小 二 乘法 拟 合 ,计算 了 震 函 
数 不 同 扩散 角度 下 的 参数 ,结果 如 表 1 所 示 。 
表 1 航 炮 冲击 波 实测 超 压 容 函 数 拟 合 结果 
Tab.1 Results of power function fitting of actual 


measured shock wave overpressure 


序号 角度/(°) Ky ag 残 差 平方 和 ( RSS) 

1 15 8 937.8 1.434 0.0006 

2 30 3755.3 1.196 0.0014 

3 45 1 444.0 0. 890 0.002 2 

4 60 3 334.4 1.142 0.0001 
= 90 4455.3 1.211 0.0006 
丁 135 8534.7 1.438 0.0063 
CO 拟 合 的 冲击 波 曲 线 与 实测 值 的 对 比如 图 6 

所 且 。 
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图 6 拟 合 曲线 与 试验 结果 的 对 比 


Fig.6 Comparison of fitting curve and test results 


2.4 超 压 时 长 


考虑 冲击 波 正 压 作 用 时 长 的 影响 ,根据 离 炮 口 
的 角度 和 距离 ,将 人 台 架 结构 的 压力 场 作用 区 进行 划 
分 。 航 炮 在 对 称 面 内 向 上 10° 角 发 射 , 按 角度 和 距离 
将 台 架 底部 蒙 皮 分 为 22 个 区 域 ,具体 分 区 与 编号 如 
图 7 所 示 。 

选择 式 (3) 函数 ,对 航 炮 炮 口 冲击 波 作用 到 台 
架 结 构 的 作用 时 长 进行 了 计算 。 以 图 7 中 的 蒙 皮 区 
域 1 为 例 , 航 炮 发 射 药 的 TNT 当量 m=0.05 kg, 炮 口 
距离 4=0.4m, 得 到 正 压 作用 时 间 1, =0.136 ms。 同 
理 ,可 以 计算 得 到 人 台 架 绽 皮 其 他 分 区 的 正 压 时 长 。 
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图 7 炮 口 冲击 波 角 度 与 压力 场 分 区 


Fig.7 Division of muzzle shock wave with 


angle and pressure field 
3 ”冲击波 仿真 


在 MCS. Patran 软件 中 ,采用 药力 试验 数据 , 建 
立 金属 材料 的 应 力 -应 变 曲线 ,如 图 8 所 示 。 
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图 8 金属 材料 实测 应 力 -应 变 曲线 
Fig.8 Measured stress-strain curves of metallic materials 

选择 Nastran 非 线性 瞬 态 SOL129 求解 器 ,对 台 
架 在 冲击 波 作用 下 的 动 响应 进行 分 析 ” 。 按 上 节 
分 区 分 别 建 立 空间 场 和 冲击 波 瞬 态 压力 载荷 ,不 同 
时 刻下 结构 表面 的 压力 分 布 变 化 如 图 9 所 示 。 

从 压力 场 的 变化 情况 可 以 看 出 , 炮 口 冲击 波 初 
始 时 刻 以 球形 形状 进行 扩散 ,最 先 到 达 离 炮 口 最 近 
的 蒙 皮 中 心 位 置 ,随后 向 外 扩散 ,同时 ,由 于 火药 燃 


气 形成 的 冠状 冲击 波 作用 , 超 压 中 心 治 发 射 方向 往 
前 扩散 ” ,由 于 正 压 时 长 的 影响 , 随 着 距离 的 增加 ， 
压力 的 持续 时 间 也 随 着 增 大 。 


三 0.03 ms 


y= TT 
ibiDib 己 一 ibm 本 me-ao 忆 一 i 


ei? 不 同时 刻 结构 表面 的 压力 场 分 布 ( 单 位 :MPa) 
CPFig.9 Pressure field distribution on structure surface 

at different moments (unit: MPa) 
设置 结构 阻尼 系数 为 0.01, 步 长 为 0.005 ms ,总 
时 半 5 ms, 对 台 架 的 动 响应 进行 了 仿真, 台 架 最 大 的 
0 10 和 图 11 所 示 。 
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图 10 台 架 最 大 变形 云图 
Fig.10 Fringes of the maximum deformation of the platform 
从 结果 可 知 , 最 大 变形 位 于 铝 蒙 皮 靠近 中 间 骨 
架 的 区 域 , 幅 值 为 9.41 mm ,对 应 的 时 刻 为 2.35 ms; 
同时 ,最 大 等 效应 力 为 162 MPa ,危险 部 位 位 于 钢 制 
蒙 皮 与 中 间 上 骨架 焊接 区 域 。 将 台 架 各 部 件 的 最 大 变 
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形 与 等 效应 力 汇总 见 表 2 所 示 。 


Patran 2014 64-Bit 08-Mar-22 15:47:55 
Fringe:cjb.Al:Time=0.00235,Nonlinear Stresses.Stress Tensor,von Mises, At Z1 
Deform:cjb.Al:Time=0.00235.Displacements,Translational 


1.62+002 
1.51+002 
1.40+002 
1.29+002 
1.19+002 
1.08+002 
9.70+001 
8.62+001 
7.54+001 
6.47+001 
5.39+001 
4.32+001 
8-002 3.24+001 

2.16+001 


1.09+001 
8.78-002 
default_Fringe: 
Max 1.62+002(@Nd 23366 


xX 
六 Min 8.78-002@Nd 18182 


六 default Deformation: 
Max 9.41+000@Nd 24442 


图 11 台 架 最 大 变形 与 最 大 等 效应 力 云 图 


Fig.11 Fringes of the maximum equivalent 


stress of the platform 
表 2 人 台 架 各 部 件 的 最 大 变形 与 最 大 等 效应 力 
Tab.2 The maximum deformation and 


equivalent stress of the platform 


部 件 底 蒙 皮 ”” 侧 蒙 皮 骨架 设备 假 件 
最 大 变形 /mm 9.41 0.6 0.38 0.39 
最 大 等 效应 力 /MPa 162 146 34.4 20.3 
峰值 时 刻 /ms 0. 235 0.251 0.235 0.27 


基于 上 述 模型 ,开展 了 参数 影响 分 析 。 选 择 台 
架 的 一 阶 垂 向 弯曲 模 态 频率 、 结 构 阻 尼 系 数 2 种 参 
数 分 别 进行 对 比 计算 ,得 到 台 架 的 最 大 变形 和 最 大 
应 力 的 变化 曲线 ,如 图 12 和 图 13 所 示 。 


台 架 最 大 虚 力 /MPa 


5 70 
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图 12 阻尼 系数 对 台 架 响应 峰值 的 影响 曲线 


Fig.12 The influence curve of damping coefficient 


on the peak response of the platform 
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图 13 台 架 固有 频率 对 台 架 响应 峰值 的 影响 曲线 


Fig. 13 The influence curve of the natural frequency of 


the platform on the peak response of the platform 
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4 结论 


1) 基于 实测 的 航 炮 冲 击 波 压力 和 经 验 公式 ,对 
炮 口 压力 场 分 布 函 数 进行 了 拟 合 ,并 考虑 正 压 的 作 
用 时 长 差异 性 ,模拟 了 真实 炮 口 冲击 波 的 传播 规律 。 
仿真 结果 与 文献 中 的 实测 规律 基本 吻合 ,验证 了 拟 
合 冲击 压力 场 的 有 效 性 和 可 靠 性 。 

2) 考 虑 金属 的 材料 非 线性 ,采用 非 线性 瞬 态 广 
法 对 某 直升机 台 架 在 航 炮 冲击 波 作用 下 的 动 响应 过 
程 进行 研究 ,结果 表明 最 大 变形 出 现在 底部 铝 蒙 皮 
的 两 侧 区 域 ,危险 部 位 位 于 中 骨架 与 钢 制 蒙 皮 的 连 
接 区 。 由 于 蒙 皮 变形 吸收 了 大 部 分 能 量 ,除了 底部 
蒙 皮 外 ,其 余 各 部 件 动 响应 较 小 ,整个 台 架 未 出 现 结 
构 破坏 。 
CD3 ) 台 架 的 响应 峰值 与 结构 阻尼 系数 成 反比 ,在 
2 和 $2350 Hz 的 范围 内 ,响应 峰值 与 台 架 的 一 阶 垂 向 栖 
则 模 态 频率 成 正比 关系 。 因 此 ,可 以 通过 减 小 对 应 
模 航 频率 或 者 增 大 结构 阻尼 来 降低 台 架 的 动 响应 。 
CD4) 通 过 上 述 方法 ,实现 了 航 炮 冲 击 波 下 的 结 检 
到 月 应 预测 ,可 以 为 相关 的 工程 结构 的 设计 提供 参 
老 肖 由 于 仿真 结果 缺少 试验 数据 对 比 ,后 续 待 开展 
相关 试验 后 对 计算 模型 进行 优化 。 
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